Millimeter Dust Emission in the GQ Lup System 
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ABSTRACT 

We present Submillimeter Array observations of the GQ Lup system at 
1.3 millimeters wavelength with 0'/4 (~60 AU) resolution. Emission is detected 
from the position of the primary star, GQ Lup A, and is marginally resolved. No 
emission is detected from the substellar companion, GQ Lup B, 0'/7 away. These 
data, together with models of the spectral energy distribution, suggest a compact 
disk around GQ Lup A with mass ~ 3 Mj up , perhaps truncated by tidal forces. 
There is no evidence for a gap or hole in the disk that might be the signature 
of an additional inner companion body capable of scattering GQ Lup B out to 
~ 100 AU separation from GQ Lup A. For GQ Lup B to have formed in situ, 
the disk must have been much more massive and extended. 

Subject headings: circumstellar matter — planetary systems: protoplanetary disks - 
stars: individual(GQ Lup) 
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Introduction 



The classical T Tauri star GQ Lup A (spectral type K7, age ~ 1 Myr) has received 
considerable attention since the detection by direct imaging in the near-infrared of a 
companion, GQ Lup B, at projected separation ~ Of! 7 (~ 100 AU at 150 pc) with a m ass 



originally claimed to be perhaps as low as the planet Jupiter (INeuhauser et al. 



2005|). 



Subsequent spectroscopy and analysis suggested a highe r mass for the companion, most 



2006 



a brown dwarf ] 


n the 


McElwain et al. 


2007 



range of 10 to 40 M 



Seifahrt et al 



2007 



lup (e.g. 



Mu grauer et al. 



Marois et al. 



2008 



2005 



Guenther 



Neuhauser et al. 



2008|), 



making it one of a small but growing class of pre-main-sequence stars with su bstellar 



companions at ~ 100 AU separations (ILuhman et al. 



2006; 



Schmidt et al. 



2009|). 



Such wide-separation, high-mass ratio companions are difficult to form in situ 
within circumstellar disks either by the standard core-accretion process o r by the 



gravi tational instability mechanism (IBoss 



2006 



Debes fc Sigurdsson 



2006 



Veras et al. 



20091 ). Observations of the disk(s) in these young systems may offer clues to their formation. 
While the GQ Lup system exhibits all of the usual pre-main-sequence disk signatures, 
including optical emission l ines due to gas accretion and thermal dust emission from infrared 



to millimeter wavelengths (INiirnberger et al. 



19971 ). it is not known whether the bulk of 



the disk mass is associated with the primary star, the substellar companion, or resides in a 
circumbinary structure. If GQ Lup B formed as a "planet" within a circumstellar disk, for 
example, then one might expect disk material to be present at large radii, beyond its orbit, 



ultimately evolvin g into a debris disk like those surro unding the recent 



orbiting HR 8799 (IMarois et al. 



20081 ) and Fomalhaut (IKalas et all 120081 ) . 



y imaged bodies 



Interferometric imaging of thermal dust continuum emission at millimeter wavelengths 
offers a way to trace the location of cool material in the GQ Lup system. In this paper, 
we present 1.3 millimeter observations with sufficient angular resolution to separate the 
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primary and the companion, which reveal a compact disk surrounding GQ Lup A. We model 
the disk emission to constrain physical properties and GQ Lup B formation scenarios. 

2. Observations 

We observed GQ Lup at 1.3 millimeters wavelength with Submillimeter Array (SMA) Q 
on 2007 June 19 as a short filler project, at the start of a set of antenna moves from the very 
extended configuration to a compact configuration. The six available antennas provided 
projected baselines from 30 to 430 meters. The two hour observation was performed in 
excellent weather conditions, with 225 GHz atmospheric opacity of 0.05, as measured at the 
nearby Caltech Submillimeter Observatory. The correlator provided 2 GHz bandwidth in 
two sidebands, with a central LO frequency of 225 GHz. Observations of the two quasars 
J1517-243 and J1454-377 were interleaved with GQ Lup (R.A. = 15 /l 49 m 12!138, DEC 
=— 35°39 / 03'/9, J2000). Passband calibration was accomplished with observations of the 
strong source 3c273. Complex gain calibration was done using J1517-243. The position 
derived for J1454-377, within 0'/2 of its catalog position, provides an empirical measure of 
the astrometric accuracy that should also hold for GQ Lup. The absolute flux scale was 
set with reference to the standard calibrator MWC349, assumed to be 1.65 Jy at 225 GHz, 
accurate to better than 20%. All of the calibration was performed with the IDL based 
MIR software. Subsequent imaging was performed using standard routines in the MIRIAD 
package. 

x The Submillimeter Array is a joint project between the Smithsonian Astrophysical Ob- 
servatory and the Academica Sinica Institute of Astronomy and Astrophysics and is funded 
by the Smithsonian Institution and the Academica Sinica. 
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3. Results and Analysis 



3.1. 1.3 mm Continuum 



Figure [T] shows the image of the GQ Lup 1.3 millimeter continuum emission, with 
synthesized beam l'/2 x 0'.'4 at position angle —28° using data from both sidebands and 
natural weighting to obtain the best signal-to-noise ratio. The beam is highly elliptical due 
to the low declination of the source combined with the limited (u, v) coverage provided by 
the short span of the observations. Nonetheless, examination of the image shows that the 
emission clearly peaks at the primary star, GQ Lup A, within the astrometric uncertainty, 
and no significant emission is detected from the substellar companion, GQ Lup B. (Note 
that any unaccounted-for orbit al motion of GQ Lup B is very small, well below the 



resolution of these observations (IMugrauer et al. 



20051 ).) The flux density determined either 



by integrating in a box around the emission in the image, or by fitting the visibilities with 
a Gaussian, is 25 ± 3 mJy (statistical error only). We identify the emission with a compact 
disk surrounding GQ Lup A. 

The SMA flux density is consistent , within the uncertaintie s, with the SEST single dish 



bolometer measurement of 38±7 mJy (INurnberger et al 



19971 ). The agreement is likely 



better than suggested by a strict comparison of the values, considering that the effective 



frequency of the ~ 50 GHz wide SES r 
the steeply rising dust spectrum (e.g. 



bandwidth is lik ely 240 to 250 GHz on account of 



Carilli et al. 



20011 ). which biases the measurement 



to a higher value (at least 20% for oc z/ 3 ). Given this consistency, we conclude that no 
significant emission has been missed due to spatial filtering by the interferometer. 

The left panel of Figure [2] shows the visibility amplitude for the SMA data as a 
function of baseline length, averaged in concentric circular annuli, appropriate for a disk 
viewed face-on. The disk inclination is unknown and cannot be determined from these data, 
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but if it is similar to th e inclination of the binary system of 27 ± 5 ( Broeg et al. 



2007 



Hiigelmeyer et al. 



20091 ). then the face-on approximation is reasonable. While the bins are 



coarse and the signal-to- noise per bin is modest (4 — 5a), the slight trend of decreasing 
amplitude at longer baselines suggests that the emission from GQ Lup is unresolved or 
marginally resolved. To estimate the atmospheric seeing effect on GQ Lup, we fit the 
calibrated visibilities of J1454-377, which has a similar angular separation from J1517-243 
as GQ Lup (albeit in a different direction), and obtain < O'.'IO fwhm. The seeing is a minor 
contributor to the visibility amplitude trend. Note that if we make the extreme assumption 
that the source emission is optically thick, then the observed flux density implies a minimum 
size. For a face-on disk with characteristic temperature 30 K, the minimum radius is 
~ 10 AU (diameter ~ 0"13). 



3.2. SED 



The right panel of Figure [2] shows the spectral energy distribution (SED) of the 
GQ Lup system, including th e new SMA datum and values from the literature. The 
optical B, V,R ,I data are from 



Herbst et al. 



(119941 ). dereddened assuming Ay =0.5 and the 



Mathisl (119901 ) e xtinction law. 



All Sky Survey ( jSkrutskie et al. 



he n ear infrared J,H,K data are from the Two Micron 



20061 ). At longer infrared wavelengths, the spectrum from 



5 to 40 u m is taken from the public archive of the Sy itzer Space Telescope l c2d' legacy 



program (jEvans et al 



2003 



Kessler-Silacci et al. 



are from the IRAS Faint Source Catalog (IMoshir et al. 



20061), and data at 12, 25, and 60 /zm 



1990) - The IRAS data appear to 



be systematically higher than the Spitzer data, more than what is readily explained by 
the presence of silicate features in the spectrum near 10 and 20 /im or by the nominal 
instrumental calibration uncertainties. 
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3.3. Disk Models 



To constrain the disk pro perties, we have calculated a series of models using the 



ra diative transfer code RADMC ( Dullemond &: Dominik 



to 



Andrews et al 



2004 ) employing an approach similar 



(120091 ) . Any constraints on the disk parameters must be considered in the 



context of the model assumptions, given the limitations of the data available. 

A natural physical model for the GQ Lup system includes a circumprimary disk with 
a sharp outer edge, as might result from tidal truncation by the secondary. Alternative 
models are possible but are not considered. In detail, we assume that the (dust) disk 
density structure is parameterized by a radial surface density power-law truncated at an 
outer radius, Rj, i.e. S = E (R/Rq)~ p , Ri n < R < Rj, and a vertical scale height pow er-law 



H = Hq(R/Rq) 1+ ^ with a "puffed-up" inner rim (as in 



Dullemond fc Dominik 



200J). The 



disk temperature structure is then calculated self-consistently, assuming stellar irradiation 



i s the only heating source. The accretion rate of 3 x 10 9 M derived by 



Hiigelmeyer et al. 



(120091 ) provides negligible heating compared to the star. For these models, we use the 
stellar spectral type and dereddened optical data to fix the stellar effective temperature, 
T* = 4060 K and luminosity, = 1.5 L , and fix R in at the dust sublimation radius of 
0.09 AU (for a sublimation temperature of 1500 K). We also fix the surface density power 
law index, p = 1. This power law ind ex is compatible with constant a irradiated accretion 



disks away from the disk boundaries^ 



D'Alcssio et al. 



observations of T Tauri star disks (Wil ner et al. 
adopt the dust opacities used by lAndrews et al 



2000 



1998), as well as resolved millimeter 



Andrews fc Williams 



20071 ). We 



(120091 ) and a standard gas-to-dust ratio 
of 100. The resulting mass opacity for (dust + gas) at 1.3 mill imeters of 0.0226 cm 2 g -1 is 



very similar to the commonly adopted value of 0.02 cm 2 g 1 of iBeckwith et al. 



(jl990f ). To 



make comparisons with the resolved millimeter data, we use a Gaussian with O'.'l fwhm 
to approximate the atmospheric seeing effect, as estimated empirically from the calibrator 
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observations. 

The three solid curves in the right panel of Figure [2] show the best-fit models to the SED 
for three values of the disk radius, Rd = 25, 50, 75 AU. Note that the silicate features near 
10 and 20 fim were excluded from the fit, as were the discrepant IRAS 12 and 25 fim data. 
The fitted values for the parameters that describe the disk vertical structure are nearly 
identical in the three models, as these parameters are constrained essentially entirely by the 
shape of the infrared SED. In each of the three models, the flaring angle is characterized by 
ip = 0.2, the scale height at 25 AU is 2.35 AU, and an inner rim scale height enhanced by 
40% is needed to account for the detailed shape of the mid-infrared spectrum relative to 
the near- infrared emission. The variations of the disk radius over the range explored have a 



negligible effect on the SED, as expected (1 Adams et al. 



19881 ). However, the corresponding 



model curves in the left panel of Figure [2] show that the disk radius has a dramatic effect on 
the 1.3 millimeter visibilities, in particular at the longer baselines that probe the relevant 
size scales. The disk models with Rd = 25 and Rd = 50 AU bracket the data points, while 
the model with Rd = 75 AU produces insufficient emission at baselines > 200 kA. 

The constraint on Rd is sensitive to the model assumptions. In particular, there is a 
well-known degeneracy between th e prescribed radial fall-off in the disk surface density 



distribution and the outer radius (jMundy et al 



19961 ). To illustrate an alternative, we 



have calculated best-fit models that assume p=1.5, a steeper surface density distribution, 
for Rd = 25, 50, 75 AU. The vertical structure parameters and SED fits in these models 
are nearly identical to the models that assume p=l. The main effect of the steeper radial 
surface density distribution is a less steep decline in the visibility function, for fixed outer 
radius, as shown by the dashed curves in the left panel of Figure [2j As expected, a slightly 
larger value of Rd is favored compared to the p=l models. A much more extreme fall-off 
of the surface density would be required to obtain Rd as large as 100 AU or larger. Given 
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the quality and noise of the millimeter data, a more complete exploration of disk model 
parameter space is not warranted. 

The disk masses in the three p=l models in Figure [2] are 0.0025, 0.0029, and 0.0032 
M q, respectively. These va lues are similar to the value of 0.0043 M Q crudely estimated 



by 



Niirnberger et al. 



of 30 K, assuming the 



(119971) from the 



Beckwith et al 



.3 mi 



limeter flux and an average dust temperature 



(119901 ) mass opacity. The uncertainty in the disk 



mass estimates is dominated by the adopted mass opacity, which depends on the dust 
properties, including grain sizes, shapes and composition, as well as the assumed interstellar 
gas-to-dust ratio. The systematic uncertainties are significant. Plausible models for dust 
properties in disks, for example including the effect of grain growth beyond millimeter sizes, 
give rise to perhaps an order-of- magnit ude range in the dust opacity at 1.3 millimeters (see 



e.g. 



D'Alessio et al 



2001 



Draine 



2006|). 



4. Discussion 



A comprehensive set of scenarios for the origi n of widely separate e 



Debes &: Sigurdssonl ( 120061 ). including 



, low mass 



companions like GQ Lup B has been considered by 
(a) formation as planet in situ at ~ 100 AU either by core accretion or by gravitational 
instability, (b) displacement of GQ Lup B to a wider orbit from a formation site much closer 
in to the primary, through stellar encounter, migration, or planet-planet scattering, and (c) 
formation as a brown dwarf by cloud fragmentation or capture. While there are arguments 
for and against each of these formation scenarios, and the true mass and nature of GQ 
Lup B remain in dispute, the observations of disk emission provide some new information. 



For in situ formation of GQ Lup B as a planet to be viable, the disk surrounding GQ 
Lup A must have extended to larger radii in the past, with sufficient mass at ~ 100 AU. 
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The detected disk is compact and low mass, only ~3 Mj up in the model with Rd = 50 AU 
described in §3.21 We can extrapolate the disk models to estimate the mass reservoir that 
might have been available for planet formation in the outer disk. For a surface density 
deceasing with radius as -R -1 , a disk extending to 250 AU, 5x larger, would have a mass 5x 
higher. In this example, the available disk mass of 12 Mj up beyond 50 AU would be only 
marginally sufficient to explain GQ Lup B, even if all of this disk mass were incorporated 
into the companion. There is no evidence for circumstellar material beyond the orbit of GQ 
Lup B, and it seems unlikely that the outer disk mass could have been substantially larger 
than this estimate (unless the adopted millimeter mass opacity were to underestimate the 
disk mass by a substantial amount). Circumstellar disks with much l arger outer radii do 



exist , e.g. IM Lup wit 



2007 



Panic et al. 



i radius ~ 900 AU in the same cloud complex (Ivan Kempen et al. 
20091 ). but these are unusual. Of course, the circumstellar disk may have 
been larger and/or more massive (and probably gravitationally unstable) at a much earlier 
evolutionary stage when the central protostar was forming. In that environment, the issue 
would be how continued infall and accretion could be quenched to allow the companion to 
remain at substellar mass. 

Mechanisms that require moving GQ Lup B outward to ~ 100 AU from closer in would 
have an effect on the inner disk, which emits primarily in the mid-infrared. The SED shows 
no significant flux deficit relative to a continuous disk that would be indicative of a gap 
or hole produced by another massive companion orbiting at < 10 AU radius that could 
have interacted with GQ Lup B to drive an outward migration or to scatter it outward. 
Moreover, such a dynamical interaction likely would leave the hypothetical inner planet 
on an eccentric orbit, disrupting the inner disk. The SED also provides no evidence for 
any major disturbance in the disk resulting from a (rare) close stellar encounter capable of 
radically changing the orbital parameters of GQ Lup B. 
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The millimeter emission from the GQ Lup system is consistent with the flux densities 



found in surve ys of classica 



1 - 100 AU flJensen et al. 



T Tauri binary systems with se parations in the range of 



1996 



Andrews fc Williams 



20051 ). with a relatively low value 



readily explained by tidal truncation of individual circumstellar disks and a paucity of 
massive cir cumbinary disks . Representative calculations of disks in binary systems by 



Artymowicz fc Lubowl (119941 ) show truncation of the circumprimary disk at < 0.5a, where 
a is the orbital semi-major axis, together with truncation of the circumsecondary disk more 
sev ere by a factor of a few . Whi le the GQ Lup mass ratio is more extreme than considered 



by I Artymowicz fc Lubowl (Il994j ) and the outer disk radius derived from the resolved 
millimeter observations is model dependent, the data are qualitatively consistent with 
expectations for tidal truncation. The SED and millimeter data show directly the presence 
of a circu mprimary disk, and the millimeter flux limit together with the suggestion of Ha 



emission ( IMarois et al. 



20071 ) could indicate a significantly smaller circumsecondary disk. 



The observations presented here limit the disk mass around GQ Lup B to be approximately 
3x less than around GQ Lup A, or < 1 Mj up . 

More sensitive dust continuum observations could show whether or not the disk around 
GQ Lup A shares a common inclination and orientation with any putative disk around GQ 
Lup B, or with the GQ Lup A-B orbit. As more substellar companions to pre-main-sequence 
stars are identified by direct imaging, deep millimeter observations will reveal more about 
their formation processes and their enigmatic relationship to planets. 



We thank A. Meredith Hughes for insightful comments, Charlie Qi for assistance with 
SMA data reduction, and Kees Dullemond for providing the radiative transfer code. We 
thank the SMA staff for scheduling the brief filler observation that provided the basis for 
this paper. 

Facility: Submillimeter Array 
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I i l i i i l i i i l i i i l 

15 h 49 m 12?3 12?2 12?1 1 2f 

RA (J2000) 



Fig. 1. — SMA image of the 1.3 millimeter continuum emission in the GQ Lup system. 
The contour levels are —2,2,4,6, ... x 2.8 mJy (the rms noise level). Negative contours are 
dashed. The star symbols mark the positions of the star GQ Lup A and substellar companion 
GQ Lup B. The millimeter emission peak is associated with GQ Lup A. The filled ellipse 

o 

indicates the l'/2 x 0'.'4, P.A. -28 synthesized beam. 
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deprojected baseline [kA] \ [/im] 



Fig. 1.— (left) The 1.3 millimeter visibility amplitude as a function of baseline length, 
averaged in three concentric circular annuli. The expected amplitude for zero signal is 
~ 5 mJy in each bin. The three solid curves are from the disk model fits to the SED 
described in §3.31 assuming p=l and disk radius values = 25,50,75 AU (in red, blue, 
and green, respectively). The three dashed curves are from the corresponding models for 
p=1.5. (right) The spectral energy distribution of GQ Lup, including the 1.3 millimeter 
SMA measurement and values from the literature. The stellar photosphere is indicated by 
the dashed line. The three solid curves are from the disk model fits, as in the left panel. 
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